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SUMMARY 


This  report  details  the  effect  of  shock- induced  conductivity  on  the 
operation  of  the  Kratz  velocity  gauge,  which  measures  the  ground  particle 
velocity  induced  by  the  shock  wave  due  to  a  large-yield  high-explosive  or 
underground  nuclear  detonation.  The  Kratz  gauge  operates  by  sensing  the 
time  rate  of  change  of  the  magnetic  field  produced  by  a  field  coil  embedded 
in  and  ideally  moving  with  the  medium.  The  shock  wave  from  a  large-yield 
detonation  can  induce  significant  conductivity  in  the  ground,  which  can  in 
turn  affect  gauge  response  both  before  and  after  the  shock  wave  passes  the 
field  coil. 

Three  models  of  increasing  fidelity  to  the  physics  of  the  advancing 
conductivity  front  are  presented  and  analyzed,  the  third  being  the  most 
complex  and  most  accurate.  Numerical  results  presented  for  model  3  indi¬ 
cate  that,  for  typical  gauge  parameter  values,  shock-induced  ground  con¬ 
ductivities  of  100  mhos/m  or  less  do  not  affect  gauge  response. 


PREFACE 


This  report  is  one  of  a  multivolume  set  comprising  the  Pacific-Sierra 
Research  Corporation  (PSR)  final  report  on  Defense  Nuclear  Agency  (DNA) 
contract  DNA001-83-C-0015 .  The  work  done  under  this  contract  spans  a 
wide  range  of  nuclear  weapon  effect  research  covering  airblast,  cratering 
and  ground  motion,  intermediate-dose  radiation,  underground  test  design 
and  development,  fire  research,  and  electromagnetic  pulse  research. 

This  volume  details  the  effect  of  shock-induced  conductivity  on  the 
operation  of  the  Kratz  particle  velocity  gauge.  The  project  officer  was 
Michael  J.  Frankel,  and  the  technical  monitor  was  David  L.  Auton. 
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SECTION  1 


INTRODUCTION 


The  Kratz  particle  velocity  gauge  [Kratz,  1976]  is  designed  to  mea¬ 
sure  the  particle  velocity  in  the  ground  induced  by  a  shock  wave  due  to 
chemical  or  nuclear  explosives.  It  has  been  successfully  used  in  several 
high-explosive  tests  as  well  as  the  nuclear  underground  tests  (UGTs) 

HURON  LANDING  and  MINI  JADE.  It  is  to  be  deployed  on  several  other  UGTs, 

including  MISTY  JADE.  The  Kratz  gauge  is  typically  used  to  measure  ve- 

3  2 

locities  ranging  from  3  *  10  m/sec  (peak)  to  8  *  10  m/sec  (peak),  decay¬ 
ing  with  a  time  constant  of  10  to  100  ysec.  The  period  over  which  mea¬ 
surements  are  made  is  a  few  time  constants. 

Both  Kratz  [1976]  and  Schlessinger  [1981]  presented  and  evaluated 
the  basic  theory  for  calibrating  the  gauge  and  quantifying  its  operation, 
but  neither  completely  modeled  the  electrical  conductivity  induced  in  the 
medium  by  the  advancing  shock  wave.  Schlessinger  demonstrated,  though 
only  approximately,  that  the  shock-induced  conductivity  would  not  signifi¬ 
cantly  influence  the  gauge  response  for  a  wide  range  of  relevant  gauge 
parameter  values.  However,  the  high  quality  of  data  recently  obtained 
from  the  gauge  [Coleman,  1982,  1983]  may  permit  us  to  identify  and  quantify 
even  the  small  changes  in  gauge  response  caused  by  the  shock-induced  con¬ 
ductivity.  A  better  theory  to  describe  those  small  signals  could  explain 
the  data  more  thoroughly. 

This  report  details  the  effect  of  shock-induced  conductivity  on 
Kratz  gauge  operation.  Because  the  problem  is  so  complex,  we  present  and 
analyze  three  models  of  increasing  fidelity  to  the  physics  of  the  moving 
conductivity  front,  the  third  being  the  most  complex  and  most  accurate. 

For  model  3,  we  give  numerical  results  that  can  be  used  to  quantify  the 
effect  of  the  shock-induced  ground  conductivity. 
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SECTION  2 


KRATZ  GAUGE  OPERATION  AND  THREE  MODELS 


Detailed  descriptions  of  the  gauge  and  its  operation  are  given  by 
Kratz  [1976]  and  Coleman  [1982].  Basically,  the  gauge  consists  of  a 
field  coil  and  several  pickup  coils.  The  field  coil  comprises  several 
turns  of  wire  embedded  in  an  aluminum  disk.  Current  is  run  through  the 
coil.  At  the  time  of  the  explosion,  the  external  current  is  turned  off 
but,  because  the  coil  is  in  aluminum,  current  persists  for  a  time  long 
compared  with  the  measurement  times  of  interest.  Thus,  the  current  in 
the  field  coil  can  be  considered  constant  in  time.  The  shock  wave  ap¬ 
proaches  and  overtakes  the  field  coil  as  shown  in  Fig.  1.  After  the 
shock  wave  passes,  the  field  coil  ideally  should  move  with  the  particle 
velocity  of  the  medium.  A  voltage  proportional  to  the  rate  of  change  of 
the  magnetic  field  and  related  to  the  field  coil  velocity  is  induced  in 
the  pickup  coils,  located  some  distance  from  the  field  coil. 


Ground 


Field  coil/current 
loop  embedded 
m  aluminum 
shell 


Shock  front  z  =  git) 


Pickup  coils 


Figure  1.  Gauge  geometry 


The  relationship  between  the  field  coil  velocity  and  the  voltage  in 
the  pickup  coil  is  modified  by  the  ambient  and  shock-induced  conductivities 
of  the  intervening  medium.  The  ambient  conductivity  allows  eddy  currents 
induced  by  the  changing  magnetic  field  of  the  field  coil  to  flow  in  the 
medium.  The  currents  produce  magnetic  field  changes  that  tend  to  cancel 
those  produced  by  the  moving  field  coil.  If  the  cancellation  is  signifi¬ 
cant,  the  signal  recorded  by  the  pickup  coil  will  be  badly  distorted. 

The  shock-induced  conductivity  can  have  two  effects.  First,  the 
pickup  coil  responds  to  the  changing  magnetic  field  produced  by  the  ad¬ 
vancing  conductivity  front  even  before  the  shock  wave  arrives  at  the  field 
coil.  Second,  the  shock-induced  conductivity  can  modify  the  gauge  re¬ 
sponse  after  the  shock  wave  passes  the  field  coil.  The  shock  wave  produces 
conductivity  between  the  field  coil  and  the  pickup  coil,  thus  modifying 
the  response  as  described  above  by  allowing  even  greater  eddy  currents  in 
the  medium. 

We  have  developed  three  models  to  describe  those  effects.  The  models 
differ  only  in  their  treatment  of  the  moving  conductivity  front;  they 
share  the  same  description  of  the  field  and  pickup  coils.  The  field  coil 
is  taken  to  be  a  current  loop  of  radius  a  centered  at  z  *  0,  r  «  0  of  a 
cylindrical  coordinate  system.  The  field  coil  carries  a  constant  current 
I  and  moves  in  the  z  direction  with  velocity  v(t)  at  time  t.  The  position 
of  the  field  coil  is  given  by 

f(t)  -  J  v(t’)  dt’  . 

— QD 

The  pickup  coil  is  a  single  loop  of  radius  R  centered  at  z  *  z^,  r  *  0 
in  the  z,  r  plane 

Model  1  assumes  that  the  coils  are  embedded  in  a  medium  of  constant 
conductivity  o  and  that  the  medium  is  moving  with  a  constant  velocity  v^ 
in  the  coordinate  system  stationary  with  respect  to  the  pickup  coil.  The 
results  of  model  1  were  presented  previously  [Schlessinger,  1981]  and  are 
only  summarized  here.  The  model  demonstrated  that  the  effects  of  the 
ambient  conductivity  were  negligible  for  all  parameter  values  of  interest 


and  that  the  shock-induced  conductivity  was  not  crucial  (for  typical 
values)  until  it  reached  about  100  mhos/m. 

Model  2  resembles  model  1  but  more  completely  describes  the  moving 
conductivity  front.  It  assumes  that  the  medium  velocity  is  a  function 
of  time — v^(t).  Both  models  assume  that  the  medium  conductivity  extends 
throughout  all  space  and  is  constant  in  space,  and  that  the  entire  medium 
is  moving  with  respect  to  the  pickup  coil.  The  results  of  model  2  are 
the  same  as  those  of  model  1  when  v^  is  held  constant.  Model  2  shows  that 
when  v^(t)  is  equal  to  v(t) — that  is,  when  the  field  coil  velocity  and 
the  medium  velocity  are  the  same — the  medium  conductivity  has  no  effect. 
Further,  when  there  are  small  differences  in  vq(0  and  v(t),  the  medium 
conductivity  begins  to  have  a  noticeable  effect,  for  typical  gauge  param¬ 
eter  values,  when  it  reaches  100  mhos/m. 

Model  3  was  devised  to  overcome  the  deficiencies  in  models  1  and  2 
and  to  be  used  for  quantifying  the  small  but  perhaps  noticeable  effects 
of  shock  conductivity.  In  model  3,  the  medium  conductivity  and  velocity 
are  allowed  to  change  as  functions  of  space  and  time.  The  medium  has 
conductivity  a  for  positions  z  <  g(t)  and  conductivity  0  for  positions 
z  >  g(t),  where  the  function  g(t)  represents  the  position  of  the  ,-:hock 
wave.  The  medium  is  moving  with  velocity  vq(0  for  z  <  g(t)  and  is  at 
rest  for  z  >  g(t).  This  model  describes  the  physics  of  the  moving  con¬ 
ductivity  front  much  more  accurately,  but  is  considerably  more  difficult 
to  describe  and  analyze  mathematically. 


SECTION  3 


ANALYSIS  OF  MODELS  1  AND  2 


-^-V- 


Models  1  and  2  are  so  similar  that  they  can  be  analyzed  together. 

The  only  difference  between  them  is  that  the  medium  velocity  vq(0  is  a 
constant  for  model  1  but  can  vary  over  time  for  model  2.  The  major  steps 
in  the  analysis  are  as  follows.  First,  we  derive  the  equations  of  the 
vector  potential  A  generated  by  the  current  density  J  in  the  field  coil. 
Second,  the  equations  for  A  are  solved  by  finding  their  Green's  function. 
Third,  we  evaluate  the  voltage  in  the  pickup  coil  generated  by  the  vector 
potential  X.  Ve  use  mks  units  throughout  this  report. 

DERIVATION  OF  BASIC  EQUATIONS 

The  Maxwell  equations  for  the  electric  field  E  and  the  magnetic 
field  B  in  a  medium  of  conductivity  cr  moving  with  velocity  vg(t)  are 

q) 


^  *  b  =  y  [  J  +  a (E  +  vQ  x  b) ]  , 


where  y  is  the  medium  permeability.  In  our  case,  the  current  density  J 
of  the  field  coil  has  a  component  only  in  the  <j>  direction  of  a  cylindri¬ 
cal  coordinate  system  and  is  given  by 


J  =  IS (r  -  a)S[z  -  f(t)J  , 


where  I  is  the  constant  current,  6(x)  is  the  Dirac  delta  function,  r  and 
z  are  the  radial  and  axial  coordinates,  respectively,  of  the  cylindrical 
coordinate  system,  a  is  the  radius  of  the  field  coil,  and  f(t)  gives  the 
position  of  the  field  coil  at  time  t. 
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To  solve  Eqs.  (1)  and  (2),  we  Introduce  the  vector  potential  A, 
satisfying 


b  =  M. 


It  follows  from  Eq.  (1)  that 


where  <t  is  the  scalar  potential.  We  choose  the  gauge  such  that 

^  •  A  +  ya4>  *»  0  . 


Because  the  current  flows  only  in  the  ip  direction,  A  has  only  a  <p  com¬ 
ponent  that  satisfies 


n2  A  3A  .  .  3A  t  , . 

v  a  -  -j  -  uo^~  yov0(t)  37  -  -  yJ  •  (4) 

r 

Here  we  have  assumed  that  the  velocity  v^Ct)  has  a  component  only  in  the 
z  direction.  Equation  (4)  must  be  solved  for  A,  which  can  then  be  used 
to  find  the  voltage  in  the  field  coil. 


GREEN'S  FUNCTION 

A  standard  method  of  solving  an  equation  like  Eq.  (4)  is  to  find 
its  Green's  function  G.  For  Eq.  (4),  that  is  the  solution  of 

^V2  -  ~  -  uo  ~  -  yovQ  ^7^  G(r ,  r',  z,  z’,  t,  t') 

*  <5(r  -  r')6(z  -  z ' )  6  ( t  -  t')  .  (5) 


Given  the  Green's  function,  the  solution  to  Eq.  (4)  is 
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A(r,  z,  t)  =  -  ]ij  dr*  dz'  dt'  G(r,  r',  z,  z  ,  t,  t')J(r',  z',  t')  . 

To  find  G,  we  express  the  delta  function  in  a  Fourier-Bessel 
representation: 


S(r  -  r')6(z  -  z’>  °  ^ ff  dk  eik(z-z,)  )  , 

where  J^(x)  is  a  first-order  Bessel  function.  We  also  express  G  in  the 
same  form  as 


G(r,  r’,  z,  z',  t,  t')  =  —  d^  eik^  *  ^ 

x  CJ1(Cr)J1(5r*)G(k,  Z,  t)  . 

The  equation  for  G  follows  from  Eqs.  (5),  (7),  and  (8): 

O  1  A  A  Jif! 

(k  +  T)G(k,  Z,  t)  +  iy0kvQG  +  yo  -  5(t  -  t')  . 

Equation  (9)  can  be  solved  to  obtain 

£<k.  5,  t)  -  -  ^  C">  exp  -  j  |  +  ikv0(t")j  dt"  |  , 

where 

0(x)  =  1  ,  x  >  0  , 

0(x)  =  0  ,  x  <  0  . 


Thus,  the  Green's  function  can  be  written  in  integral  form  using  Eqs. 
(8)  and  (10),  which  can  be  further  reduced  by  performing  some  of  the 


integrations.  The  £  integration  can  be  performed  using  the  result 
[Gradshteyn  and  Ryzhik,  1980] 


where  I^(x)  is  the  Bessel  function  of  the  second  kind  and  first  order. 
Also,  the  k  integral  can  be  done  using  the  result 
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EVALUATION  OF  VOLTAGE  IN  PICKUP  COIL 

The  voltage  V  in  the  pickup  coil  is  determined  by  the  rate  of  change 
of  the  magnetic  flux  linking  the  coil.  In  terms  of  the  vector  potential 
A,  the  voltage  is  given  by 


V(t)  =  2ttR  A(R,  zQ,  t) 


r 

■  ■ 
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w 
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Using  Eq.  (6),  this  expression  becomes 


V(t)  -  2ttR 


TtffJ-  yG(R*  zo*  c’  r’’  z’’  c,) 


x  J(r\  z',  t')  dr'  dz’  dt'  . 


It  follows  from  Eqs.  (3)  and  (11)  that  V(t)  is 


v(t) .  wtf  £  /*»• 


r. _ ms _ / 

L  ^(t  -  t’)  \ 


(z  -  f(t’)  -  D]2  +  (R  -  a)2 


Taking  the  indicated  derivatives  and  changing  variables  results  in  the 
final  expression  for  the  voltage: 
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V(t)  *  ^  Ray I  y  dx  /x  [zQ  -  f(y)  -  D][v(y)  +  vQ( t)  -  vq(Y)] 
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SECTION  4 


ANALYSIS  OF  MODEL  3 


Model  3  encompasses  models  1  and  2  but  also  includes  the  effects  of 
the  moving  shock-induced  conductivity  front.  In  model  3,  the  field  coil 
is  located  at  position  z  *  f(t),  as  before,  and  the  shock  front  at 
z  *  g(t).  The  medium  has  conductivity  a  and  velocity  v^(t)  for  positions 
behind  the  shock  wave  [z  <  g(t)];  the  conductivity  and  velocity  are  zero 
for  positions  in  front  of  the  shock  wave  [z  >  g(t)].  We  label  as  region  1 
the  shocked  region,  where  a  ^  0;  and  as  region  2  the  unshocked  region, 
where  a  *  0.  The  fields  that  exist  before  the  interface  (shock  wave)  has 
passed  the  field  coil  are  labeled  with  a  superscript  b. 

Before  the  interface  reaches  the  field  coil,  the  field  in  region  1 
is  a  solution  to 


(v2  -  7  - ^  -  wavo  £)  a5  =  0 


This  solution  can  be  represented  as 


A1  °  "  yI0  /-I  dk2 


x  exp 


+  ik2D  -  ik2Z  -  - 


(kl  +  k2) 


ua 


aj(klf  k2) 


In  region  2,  before  the  interface  reaches  the  field  coil,  the  field 
is 


.  h  b  .  b 
A2  =  A0  +  a2  • 


where  A 


b 

n 


is  the  field  produced  by  the  field  coil. 


(15) 
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(16: 


and  a. 


Is  a  solution  of 


(-•  ■  ?) 


ab  ■  0 


This  expression  can  be  represented  as 


=  -  yIQ  j dk  e~kz  J^krjal^k) 


The  Green's  function  G  is  given  by 


Gb(r,  r',  z,  z')  =  -|y*dC  exp  (-  £|z  -  z’|)  J^CO  J^Cr  ’ )  .  (17 


The  functions  ab  and  ab  are  determined  by  requiring  that  A  and  3A/3z 


be  continuous  at  z  *  g(t) — that  is,  across  the  shock  front.  Imposing 
those  conditions  yields  the  following  equations  for  a^,  a2: 
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dk^  exp 


+  ik2D  -  ik2g  - 


(kl  +  k2  ) 


yo 
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-  a  exp  [-  k1(f  -  g)]  J1(k1a) 
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and 
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dk2  exp 
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2  exp  (-  k.g)  a2(k  )  . 
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with  G  given  by  Eq.  (11),  and 


UIo  /dki 


dk2  J^(k^r) 


x  exp 


ik2D 
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VO 
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The  field  in  region  2,  after  the  interface  has  passed  the  field  coil, 
is  given  by 


=  -  m 


o  / 


dk  e 


-kz 


J1(kr)a2(k) 


(20) 


The  quantities  a^  and  a2  are  determined,  as  before,  by  requiring  A  and 
3A/3z  to  be  continuous  at  z  *  g(t).  Imposing  those  conditions  yields 
the  following  equations  for  a^,  a2: 


17 
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Equation  (21)  can  be  solved  for  a^(k^,  k2>  and  used  in  Eq .  (22)  to  ob¬ 
tain  a2(k^,  t). 

The  voltage  in  the  pickup  coil  is  then  obtained  using  Eqs.  (12)  and 
(15),  (16),  (17),  or  (20). 

Numerical  solution  of  Eqs.  (18),  (19),  (21),  and  (22)  is  difficult 
because  a^  is  a  sharply  peaked  function  having  a  maximum  af  «  -  ik2> 
Thus,  an  appropriate  method  of  solution  is  used  here.  Equation  (18)  can 
be  satisfied  only  if 


-  kjg  -  ik2(g  -  D)  - 


(kl  *  k2  ) 


^ 1  =  MO(g  -  P)  _  ( [\  +  ua(g  -  D)12  +  1*0D  l 
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For  small  values  of  cr,  this  expression  is 


2  ,  uaD 

kx  "  1  "  2kxt  • 


For  large  values  of  a,  it  is 


-r-=-  « - - 

\  8  ~  D 


Assuming  a^(k^,  k2^  ~  al*^k2  ~  k^) ,  then  Eqs.  (18)  and  (19)  can  be  solved 
to  yield 


a  Rul  yaD _ 1 _ 
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for  0  small  and 


2„ 

a  RylD 


2  4(2g  -  D)  (z  +  £  _  2g)3 


(23b) 


for  0  large.  For  our  problem,  a  is  almost  always  small  (<10J  mhos/ra) . 
Thus,  using  Eq.  (23a),  we  evaluate  the  voltage  in  the  pickup  coil  before 
the  field  coil  starts  to  move  (i.e.,  g  <  f)  as 


For  the  case  when  the  shock  wave  has  passed  the  field  coil,  the  con 
tribution  of  a^  to  Eq.  (22)  is  negligible.  Thus, 


After  the  small-argument  approximation  for  the  Bessel  functions  has  been 
used,  one  integration  can  be  performed  and  the  result  put  into  a  more 
convenient  form: 
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SECTION  5 
NUMERICAL  RESULTS 
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To  obtain  numerical  results,  we  evaluate  the  expressions  in  Eqs.  (14), 
(24),  and  (26).  The  velocities  were  parameterized  as  zero  for  t  <  0  and  as 


'0(t)  ■  vo  [ew  (-  - exp  (-  ^ )]  > 


-  4xp  [- ^r1]  - exp  [-  iV2i]}- 


v(t) 


v  (c)  *  v 
s  s 


[“?  (•  t)  -  e*p  (•  t)] 


for  t  >  0  and  all  velocities  are  zero  for  their  arguments  less  than  zero. 
The  distances  are  taken  as  zero  for  t  <  0  and  as 


°(t)  ■  v0  |[l  -  “P  (-  tj\  T1  ■  [L  -  e*P  (-  ^)]  T2  |  ' 


f(t)  =  f0  +  V 


I3i  - [-  -  i1  •  “p[-  • 


:<0.vsj[l-e»p(-^)]x4-[l-exp(-i)] 


•  4 


for  t  >  0. 


We  evaluated  Eqs.  (14),  (24),  and  (26)  using  the  following  values 

for  f0’  z0’  V’  V  V  T1  throu8h  t4»  co*  and  a: 

f0  =  8^t0^  =  0-25  m* 

Zq  =  0.75  m, 

—  3 

v  =  10  ra/sec, 

vQ  =  v  exp  =  1*63  x  102  m/sec, 

—  3 

v  =  3  x  10  m/sec, 
s 

=  2  X  10  4  sec, 

=  10  ^  sec, 

.  -3 

=  10  sec, 

tQ  =  9.8  x  io-5  sec, 
a  =  0,  10,  100  mhos/ra. 

These  parameter  values  are  similar  to  those  appropriate  for  MINI  JADE 
and  MISTY  JADE.  The  results  fcr  these  values  are  plotted  for  models  2 
and  3  in  Figs.  2  and  3,  respectively. 

Knowing  the  value  of  the  shock-induced  conductivity  and  using  the 
theory  developed  above,  it  is  possible  to  determine  the  particle  velocity 
from  the  gauge  response.  However,  the  value  of  the  shock-induced  conduc¬ 
tivity  is  rarely  well  known.  Thus,  the  gauge  should  be  used  to  measure 
particle  velocity  only  in  situations  where  the  effect  of  the  conductivity 
does  not  depend  strongly  on  the  value  of  the  conductivity. 

Figures  2  and  3  show  that,  for  typical  gauge  parameter  values,  shock- 
induced  conductivity  just  begins  to  affect  the  gauge  response  at  100  mhos/m 
We  thus  conclude  that,  for  measurements  of  MINI  JADE  and  MISTY  JADE  phe¬ 
nomena,  shock-induced  conductivity  has  little  effect  on  gauge  response  as 
long  as  it  is  less  than  100  mhos/m. 

Figure  3  shows  the  magnitude  and  time  dependence  of  the  precursor 
signal  caused  by  the  advancing  conductivity  front  before  the  field  coil 
begins  to  move.  It  has  been  suggested  that  the  precursor  signal  could  be 

f 


Relative  gauge  output 


used  Co  deduce  the  medium  velocity  before  the  field  coil  begins  to  move. 
However,  Fig.  3  indicates  that,  when  the  gauge  is  designed  to  measure  the 
signal  from  the  moving  field  coil  (and  consequently  O  <  100  mhos/m),  the 
precursor  signal  is  probably  too  small  to  be  reliably  used  to  measure  the 
medium  velocity. 

If  one  is  willing  to  design  the  gauge  to  measure  only  the  precursor 
signal,  then,  according  to  Eqs.  (24a)  and  (24b)  for  the  gauge  response, 
in  certain  ranges  of  conductivity,  the  gauge  could  measure  either  the 
shock-induced  conductivity  or  the  medium  velocity.  That  is,  Eq.  (24a) 
shows  that,  for  0  <  1000  mhos/m,  the  gauge  response  is  linear  in  the  con¬ 
ductivity,  so  precursor  signal  measurements  along  with  independent  mea¬ 
surements  of  the  medium  velocity  and  shock  velocity  could  be  used  to 
determine  the  medium  conductivity.  Further,  Eq.  (24b)  shows  that,  for 
a  >  1000  mhos/m,  the  precursor  signal  is  once  again  independent  of  the 
shock-induced  conductivity.  Thus,  for  very  large  conductivities  (>1000 
mhos/m),  the  precursor  signal  could  be  used  to  measure  the  medium  velocity. 

We  conclude  that  the  Kratz  particle  velocity  gauge  is  not  sensitive 
to  shock-induced  conductivity  for  conductivities  less  than  about  100  mhos/m. 
For  gauges  designed  to  measure  the  velocity  of  the  field  coil,  the  pre¬ 
cursor  signal  is  probably  too  small  to  be  useful.  However,  gauges  can  be 
designed  to  take  advantage  of  the  precursor  signal  and  measure  either 
the  shock-induced  conductivity  (for  O  <  1000  mhos/m)  or  the  medium  ve¬ 
locity  (for  a  >  1000  mhos/m). 
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APPENDIX 


LIST  OF  SYMBOLS 


a  »  radius  of  field  coil 

a^  *  vector  potential  in  region  2  (which  has  no  shock-induced 
conductivity)  after  shock  wave  has  passed  field  coil 

contribution  to  vector  potential  in  region  2  (which  has 
no  shock-induced  conductivity)  caused  by  presence  of 
conductivity  front  before  shock  wave  has  passed  field 
coil,  i.e. ,  A&  =  Aj>  +  a& 

_  b  b 

Bessel  transforms  of  a^,  a^,  a^,  a^ 

vector  potential 

Aq  *  vector  potential  of  field  coil  in  homogeneous  moving 
medium  with  shock-induced  conductivity 

A^  *  contribution  to  vector  potential  in  region  1  (which  has 
shock-induced  conductivity)  caused  by  presence  of  con¬ 
ductivity  front  after  shock  wave  has  passed  field  coil 

Aq  =  vector  potential  of  field  coil  in  homogeneous  medium  of 
zero  conductivity 

A^  »  vector  potential  in  region  1  (which  has  shock-induced 

conductivity)  before  shock  wave  has  passed  field  coil 

A^  =  vector  potential  in  region  2  (which  has  no  shock-induced 
conductivity)  before  shock  wave  has  passed  field  coil 

b  *  superscript  indicating  fields  that  exist  before  shock 
wave  has  passed  field  coil 

B  *  magnetic  field 

D  =*  time  integral  of  medium  velocity 

E  *  electric  field 
f(t)  *  axial  position  of  field  coil 
g(t)  *  axial  position  of  shock  wave 

G  *  Green's  function 

G^  s  Green's  function  for  field  coil  in  medium  of  zero 
conductivity 


I  =  current  in  field  coil  (assumed  constant) 


j  =  current  density 

r  =  radial  coordinate  of  cylindrical  coordinate  system 
R  =  radius  of  pickup  coil 

U(a,  b,  x)  =  confluent  hypergeometric  function  of  second  kind 
(singular  at  x  ■  0) 

v(t)  *  field  coil  velocity 

Vg(t)  “  medium  velocity 

vg  *  axial  velocity  of  shock  wave 

V  =•  voltage  in  pickup  coil 

z  *  axial  coordinate  of  cylindrical  coordinate  system 
Zq  *  axial  position  of  pickup  coil 
6(x)  ■  Dirac  delta  function 

$  *  azimuthal  angle  of  cylindrical  coordinate  system 
♦  ■  scalar  potential 

tp(af  b,  x)  *  confluent  hypergeometric  function  of  first  kind  (regu¬ 
lar  at  x  *  0) 

O  *  medium  conductivity 

y  *  medium  permeability  (taken  to  be  that  of  free  space) 
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